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1.0  INTRODUCTION 


The  purpose  of  this  document  is  to  prescribe  the  minimum  requirements  for 
measuring  environment?.!  parameters  that  affect  the  evaluation  of  proposed  radai  systems. 
These  environmental  factors  must  be  properly  measured  to  ensure  fair  appraisals  of 
competing  designs  and  to  assess  expected  long-term  performance.  The  procedures  for 
short-term  in-situ  measurements  pertinent  to  microwave  propagation  and  sea  clutter  are 
addressed  in  sections  2.0  through  4.0.  Measurement  accuracies  are  specified  and  several 
types  of  instruments  are  discussed  in  these  sections.  Specific  manufacturers  and  retailers  are 
mentioned  as  sources  for  various  types  of  instruments  and  equipment  in  section  5.0.  This 
listing  is  not  all-inclusive  and  is  not  intended  to  be  an  endorsement  of  those  companies  or 
products.  Section  6.0  is  a  glossary  of  terms  used  in  instrument  specifications.  Appendixes  A 
through  F  are  included  as  background  material  for  determining  refractive  structure  from 
surface  meteorological  measurements,  f  >'  p>  j  — - 

Environmental  measurements  relevant  to  clutter  from  clouds,  rain,  land,  birds, 
and  insects  are  not  addressed.  Environmental  measurement  specifications  relevant  to  the 
performance  of  infrared  systems  are  set  forth  in  a  separate  NATO  Anti-Air  Warfare 
System  (NAAWS)  document. 

Previous  work  performed  under  the  NAAWS  project  has  shown  that  certain  envir¬ 
onmental  factors  strongly  affect  overall  radar  system  performance.  For  shipborne  radars 
operating  at  frequencies  above  L-band,  the  evaporation  duct  is  the  dominant  propagation 
mechanism  affecting  the  maximum  detection  ram  ^  horizon-search  radars.1  In  addition, 
surface-based  ducts  may  also  affect  the  detection  c.  •'•a*  ’ities  of  horizon-search  radars. 
Ducting  can  also  increase  sea  clutter  return  within  a i.  oeyond  the  normal  horizon,  and 
surface-based  ducts  can  enhance  land  clutter  return  from  greatly  extended  ranges.2  Sea 
clutter  return  is  also  dependent  upon  sea  state,  which  in  turn  is  dependent  upon  wind 
speed,  wind  direction,  fetch,  and  duration. 
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2.0  MARINE  ATMOSPHERIC  BOUNDARY  LAYER 


The  marine  atmospheric  boundary  layer  is  the  coupling  region  between  the  free 
atmosphere  and  the  ocean.  Air-sea  interaction  effects  occur  at  the  bottom  of  the  boundary 
layer,  within  the  surface  layer,  with  the  transport  of  momentum,  heat,  and  moisture.  The 
surface  layer  is  a  few  tens  of  meters  thick.  Above  the  surface  layer,  buoyant  and  mechanical 
turbulence  tend  to  yield  a  well-mixed  layer  within  which  gradients  of  meteorological 
parameters  are  constant.  The  thickness  of  the  mixed  layer  varies  from  a  few  hundred  to  a 
thousand  meters.  Often,  the  mixed  layer  is  capped  by  an  interfacial  layer  that  provides  a 
transition  between  the  mixed  layer  and  the  free  atmosphere.  The  interfacial  layer  may  be 
quite  weak  or  even  nonexistent.  If  it  is  present,  this  layer  is  characterized  by  an  increase  in 
temperature  and  a  decrease  in  moisture.  Idealized  profiles  of  temperature  and  moisture  in 
the  boundary  layer  are  shown  in  Fig.  1,  which  is  adapted  from  Davidson  et  al.3 


Figure  1.  Idealized  profiles  of  potential  temperature,  6,  and  specific 
humidity,  q,  in  a  well-mixed  boundary  layer.  (Adapted  from  Ref.3.) 
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3.0  UPPER  AIR  MEASUREMENTS 


Measurements  of  the  profiles  of  Fig.  1  in  the  mixed  layer  and  above  are  usually 
done  with  a  radiosonde  carried  aloft  by  a  free  balloon.  The  radiosonde  measures  pressure, 
temperature,  and  (usually)  relative  humidity  as  it  ascends  through  the  atmosphere  and 
telemeters  these  data  to  a  ground-station  receiver.  The  output  of  the  ground  station  may  be 
either  digital  or  analog.  The  analog  data  must  be  converted  to  engineering  units  of  pressure 
( P  in  millibars),  temperature  (Tin  Celsius),  and  relative  humidity  ( RH  in  percent).  These 
data  can  then  be  converted  to  refractivity,  N 


77.6  P  .  e 

N=  -  +  3.73  x  I05  — r 

T  T 2 


or  modified  refractivity,  M 


M-  A  +  0.157A 


where  P  is  pressure  in  millibars,  T  is  temperature  in  Kelvin,  e  is  water  vapor  pressure  in 
millibars,  and  h  is  height  in  meters.  The  ambient  water  vapor  pressure,  e ,  is 


RH 


where  RH  is  relative  humidity  in  percent  and  saturation  water  vapor  pressure,  es,  in  milli¬ 
bars,  can  be  determined  from 


=  6. 105  exp  25.22 


T-  273.2 


5.31  In 


where  Tis  air  temperature  in  Kelvin  and  In  is  the  natural  logarithm. 

3.1  REFRACTIVE  STRUCTURE 

Refractive  gradients  have  been  qualitatively  characterized  by  their  effects.  The 
standard  atmosphere  gradient  determined  from  climatological  data  has  been  defined  as 
1 18  Mj  km.  This  gradient  causes  an  electromagnetic  wave  to  bend  slightly  toward  the 
surface  of  the  earth  so  that  the  radar  horizon  is  slightly  greater  than  the  geometric  horizon. 
Of  course,  this  value  varies  with  meteorological  conditions;  gradients  different  from 
1 18  Ml  km  that  do  not  have  dramatic  effects  on  propagation  are  referred  to  as  normal. 
Gradients  less  than  normal  cause  electromagnetic  waves  to  bend  even  more  toward  the 
earth’s  surface  and  are  referred  to  as  superrefractive.  The  extreme  case  of  superrefraction  is 
called  trapping,  when  electromagnetic  waves  are  bent  downward  with  a  curvature  greater 
than  that  of  the  earth.  The  trapping  layer  thus  forms  a  duct  that  can  behave  very  much  like 
a  waveguide.  Gradients  much  greater  than  118  Ml  km  are  called  subrefractive  and  result  in 
electromagnetic  waves  being  bent  away  from  the  earth’s  surface.  Table  1  summarizes  the 
refractive  gradients. 
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Tabic  I.  Relation  of  A-  and  A/-gradients. 


Gradient  Type 

A-Gradient 

A/-Gradient 

Trapping 

<157  A/ km 

<0  M/km 

Superrefractive 

-  157  to  -  79  A/km 

0  to  79  Ml  km 

Normal 

-79  to  0  A/  km 

79  to  157  Mj  km 

Subrefractive 

>  0  A/  km 

>157  Ml  km 

The  temperature  and  moisture  gradients  in  the  interfacial  layer  usually  result  in  a 
superrefractive  layer.  If  the  gradients  in  the  interfacial  layer  are  sufficiently  strong,  a  trap¬ 
ping  layer  will  be  formed.  If  the  trapping  layer  occurs  well  above  the  surface,  an  elevated 
duct  will  be  formed  (Fig.  2a).  If  the  gradients  in  the  interfacial  layer  are  stronger,  or  the 
trapping  layer  occurs  at  a  lower  height,  a  surface-based  duct  will  be  formed  (Fig.  2b).  These 
two  types  of  duct  are  formed  under  stable  atmospheric  conditions  and  vary  slowly  in  range 
and  time.  Under  unusual  conditions,  the  interfacial  layer  may  occur  at  the  surface,  creating 
a  surface-based  duct  from  a  surface  trapping  layer  (Fig.  2c)  This  situation  typically  occurs 
downwind  of  a  large  land  mass  and  is  often  characterized  by  a  strong  variation  of  the 
refractive  structure  with  range  and  time. 

An  elevated  duct,  whose  base  is  more  than  a  few  hundred  meters  above  the  surface, 
will  have  little  impact  on  surface-to-surface  propagation.  A  surface-based  duct  from  an 
elevated  refractive  layer  will  have  little  impact  within  the  normal  radar  horizon  but  may 
provide  detection  capability  and  significant  clutter  return  from  well  beyond  the  horizon. 
The  surface-based  duct  from  a  surface  trapping  layer  will  be  the  dominant  propagation 
mechanism  when  this  condition  exists.  Radar  propagation  and  clutter  signals  will  be 
dependent  upon  the  gradient  and  thickness  of  the  trapping  layer. 

Although  trapping  layers  and  the  associated  ducts  give  rise  to  the  most  dramatic 
effects,  superrefractive  and  subrefractive  layers  may  also  affect  propagation  and  clutter. 
Craig  et  al.  discuss  these  various  meteorological  effects  in  detail.4 
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ALTITUDE 


(a)  Elevated  duct.  (b)  Surface-based  duct  from 

an  elevated  trapping  layer. 


(c)  Surface-based  duct  from 
a  surface  trapping  layer. 


Figure  2.  Idealized  modified  refractivity  profiles. 
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3.2  MEASUREMENTS  AND  PROCEDURES 


When  the  vertical  profile  of  refractivity  is  being  measured,  the  absolute  accuracy  of 
the  data  is  not  as  important  as  an  accurate  measurement  of  the  gradients.  Sensor  lag  tends 
to  reduce  gradients  and  increase  the  altitudes  at  which  the  changes  in  gradient  appear.  The 
lag  errors  are  not  too  serious  at  normal  ascension  rates  but  can  be  very  significant  in  the 
iirst  several  tens  of  meters  above  the  surface  if  the  sensor  has  been  abnormally  heated  or 
cooled  just  before  launch/  This  mainly  affects  the  radiosondes  with  carbon  hygristor 
humidity  elements,  which  have  a  relatively  large  thermal  mass,  but  other  types  of  humidity 
sensors  and  temperature  sensors  can  also  be  affected  by  surface  or  platform  heating.  Care¬ 
ful  preparation  and  handling  of  the  instrument  package  before  launch  can  minimize  these 
problems. 

Table  2  shows  the  minimum  specifications  for  radiosonde  sensors.  For  conditions 
typical  of  the  NAAWS  area  of  interest,  these  sensor  specifications  will  nominally  provide 
1  to  2  A-units  resolution  as  calculated  from  the  root  of  the  sum  of  the  squares: 

A N  =  [  (aAP)2  +  (6A7)2  +  (c&RH)2  ]V: 
where  a  =  dN/dP,  b  =  dN/dT,  and  c  =  dN/dRH.b 

Table  2.  Minimum  specifications  for  radiosonde  sensors. 


Parameter 

Precision 

Resolution 

Range 

Temperature 

0.5°C 

0.1°C 

- 10  to  40°  C 

Relative  humidity 

5% 

i  % 

20  to  100% 

Pressure 

3  mb 

0.5  mb 

850  to  1050  mb 

T  he  two  main  types  of  instruments  that  can  be  used  to  obtain  refractivity  profiles  of 
the  boundary  iayer  are  the  radiosonde  and  the  refractometer.  Environmental  measurements 
using  the  refractometer.  while  practical,  are  neither  simple  nor  inexpensive.  It  used,  the 
refractometer  should  have  a  precision  of  ±2  A-units. 

Radiosondes  are  commonly  borne  aloft  by  balloons.  Radiosondes  packaged  to  be 
launched  from  an  aircraft  (dropsondes)  are  also  available.  In  recent  years,  the  technology 
for  launching  a  radiosonde  package  with  a  small  rocket  has  been  developed;  the  rocket 
boosts  the  package  to  a  given  altitude,  where  it  is  ejected  to  descend  by  parachute  as  it 
makes  measurements. 

Balloon-borne  radiosondes  are  relatively  inexpensive  at  a  cost  of  approximately 
S100  (U.S.)  per  launch  for  sonde,  balloon,  and  gas.  Ground-station  equipment  is  consider¬ 
ably  more  expensive,  depending  upon  the  features  of  the  equipment.  This  equipment  may 
be  available  for  lease.  Radiosondes  carried  aloft  by  an  aerostat,  under  the  control  of  a 
winch,  are  also  an  option.  Detailed  sounding  procedures  will  be  specified  in  manufacturer’s 
instructions.  Procedures  applicable  to  NAAWS  include  selecting  an  appropriate  site  (or 
sites)  for  launch  and  handling  that  minimizes  spurious  sensor  readings.  Spurious  sensor 
readings  can  arise  from  radiative  and  conductive  heating  from  surrounding  structures  or 
from  failing  to  acclimatize  the  radiosonde  to  ambient  conditions  after  removing  it  from  a 
heated  or  air-conditioned  space. 
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Aboard  a  boat  or  ship  the  radiosonde  should  be  launched  from  as  low  an  altitude 
as  possible  on  the  windward  side  of  the  bow  or  stern  in  order  to  minimize  the  heat  island 
effect  of  the  ship.  Care  should  be  taken  to  acclimatize  the  sonde  to  ambient  conditions 
before  releasing  it;  that  is.  the  sensors  should  be  given  at  least  15  seconds  to  adjust  to  the 
outside  conditions  before  launch.  Similar  procedures  are  appropriate  for  launch  sites  on 
coastal  structures. 

If  the  radiosonde  is  launched  from  shore,  the  launch  site  should  be  as  close  to  the 
water  as  possible,  and  precautions  should  be  taken  to  minimize  radiative  and  conductive 
heating  effects  from  parking  lots,  buildings,  and  other  structures.  Near  active  airways, 
coordination  with  air  control  personnel  may  be  required. 

If  the  wind  at  the  time  of  the  radiosonde  observation  is  such  that  the  air  has  had 
a  long  overwater  trajectory  (>100  km),  the  boundary  layer  should  be  representative  of 
open-ocean  conditions.  The  assumption  of  horizontal  homogeneity  usually  will  be  valid, 
and  one  sounding  will  generally  be  representative  of  conditions  over  the  path  of  interest  to 
NAAWS.  However,  if  the  wind  is  offshore  or  alongshore,  propagation  conditions  over  the 
propagation  path  are  likely  to  be  inhomogeneous,  and  a  minimum  of  two  soundings  should 
be  taken  along  the  path  within  one-half  hour  of  each  other. 

Balloon  ascent  rates  should  be  in  the  range  of  200  to  300  m  min.  Slower  ascent 
rates  yield  a  higher  vertical  resolution  in  the  data  but  may  not  provide  adequate  ventilation 
of  the  sensors.  Faster  ascent  rates  provide  good  sensor  ventilation,  but  some  vertical  resolu¬ 
tion  is  lost. 

All  data  collected  during  the  sounding,  whether  digital  or  analog,  should  be 
retained.  Sounding  systems  that  provide  automatic  data  analysis  use  criteria  for  selecting 
significant  data  that  are  oriented  toward  fulfilling  meteorological  data  requirements;  these 
criteria  may  be  too  lenient  for  detailed  propagation  assessment. 
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4.0  SURFACE  MEASUREMENTS 


4.1  SURFACE  LAYER 

The  atmospheric  surface  layer  is  usually  dominated  by  mechanical  and  buoyant 
turbulence.  The  average  profile  of  wind,  temperature,  humidity,  or  any  conservative  scalar 
is  determined  by  the  turbulent  motions.  Roll  presents  an  overview  of  the  wind,  tempera¬ 
ture.  and  moisture  fields  in  the  first  few  meters  above  the  sea  surface.7  Direct  measurements 
of  the  profiles  in  the  surface  layer  over  the  ocean  are  quite  difficult  and  essentially  imprac¬ 
tical  on  an  operational  basis.  To  circumvent  this  difficulty,  semiempirical  relationships 
between  the  profiles  and  the  fluxes  of  momentum,  heat,  and  moisture  have  been  developed. 

4.1.1  Refractive  Structure 

In  meteorology,  it  is  convenient  to  work  with  parameters  which  are  invariant  with 
height  in  homogeneous  air.4  These  are  often  referred  to  as  potential  parameters  and 
conservative  properties  with  respect  to  adiabatic  processes.  For  radio  propagation,  the 
profile  of  potential  refractivity  is  of  interest,  with  potential  refractivity  defined  by 


77.6  P„ 

ft 


+  3.73  » 


I05 


where  Pn  -  1000  mb.  ft  is  potential  temperature  in  Kelvin,  and  ep  is  potential  vapor  pres¬ 
sure  in  millibars.  Further 


ft  -  T 


0 


and 


where  T  is  temperature  in  Kelvin.  P  is  pressure  in  millibars,  and  e  is  water  vapor  pressure 
in  millibars.  Over  the  ocean,  the  profile  of  potential  refractivity  is  usually  such  that  trap¬ 
ping  gradients  exist  in  the  first  few  meters  of  the  surface  layer,  and  a  shallow  duct,  called 
the  evaporation  duct  because  of  the  predominant  influence  of  evaporation  in  causing  the 
trapping  gradient,  is  found.  The  parameter  that  best  characterizes  propagation  is  the  evap¬ 
oration  duct  height,  defined  as  the  height  at  which  the  gradient  of  potential  refractivity 
reaches  the  critical  value,  calculated  by  rav  theory,  which  will  cause  a  ray  launched  horiz¬ 
ontally  to  have  a  curvature  exactly  equal  to  the  curvature  of  the  earth.  In  terms  of  modified 
refractivity.  M.  the  cvapoiation  duct  height  is  the  height  in  the  surface  layer  at  which  the  M 
profile  reaches  its  minimum  value  and  the  .^/-gradient  is  zero.  Figure  3  shows  an  example 
of  an  XI  profile. 

There  are  several  semiempirical  techniques  to  relate  surface  profiles  to  meteorolog 
ical  measurements  that  can  be  simply  made  at  sea.  Jeske  examined  several  schemes  and 
correlated  evaporation  duct  height  to  propagation  measurements.8  10  His  work  is  the  basis 
for  the  duct  height  formulation  that  follows;  this  formulation  is  developed  in  more  detail  in 
Appendixes  A.  B.  and  C.  A  FORTRAN  77  program  and  test  case  calculations  are  provided 
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285  M  305 

MODIFIED  REFRACTMTY 

Figure  3.  Example  of  a  modified  refractivity 
profile  for  a  10-m  evaporation  duct.  At  the 
surface,  M  =  300;  z  is  altitude. 

in  Appendixes  D  and  E,  respectively.  Appendix  F  develops  the  equations  for  the  M-pr ofile 
in  terms  of  duct  height  and  stability  for  use  in  propagation  models.  Rotheram  has  devel¬ 
oped  a  similar  formulation  in  a  slightly  different  manner." 

The  following  six  steps  specify  the  calculation  of  the  evaporation  duct  height: 

1.  The  four  bulk  meteorological  measurements  required  are  air  temperature,  Ta , 
and  sea  temperature,  Ts,  in  Celsius;  relative  humidity,  RH ,  in  percent;  and  wind  speed,  u, 
in  knots.  Temperatures  are  converted  to  Kelvin  by  Tak  =  Ta  +  273.2  and  Tsk  =  Ts  +  273.2.  If 
the  wind  speed  is  less  than  0.01  knot,  the  evaporation  duct  height  is  set  to  0  meters. 

2.  The  bulk  Richardson  number  is  calculated  by 
/?/,  =  369  z,  ( Ta-Ts )/  (Tak  u2) 

where  z,  is  the  reference  altitude  in  meters.  Rih  is  restricted  to  be  no  greater  than  1. 

3.  Assuming  that,  in  the  surface  layer,  Tak  0,  e  «  e  ,  and  neglecting  the  effect 
of  salinity  on  vapor  pressure,  calculate  potential  refractivity  for  the  air.  NA,  at  height  z, 
and  for  the  air  immediately  in  contact  with  the  sea  surface,  N 


77.6 

f 

4810 

A  -r  ( 

1  ak 

1000  + 

Tak 

77.6 

/ 

4810 

7  —  | 

[1000  + 

c  1 

Tik 
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where  the  ambient  vapor  pressure  of  the  air,  e,  is  calculated  by 


with 


RH 

e  -  -  e, 

100 


T  -  27”}  2  T 

es  -  6.105  exp  (  25.22  ---  _ — —  -  5.31  In  ak 


Tok  273.2 

and  the  vapor  pressure  at  the  sea  surface,  ^0,  is 


T  -  27”}  2  T 

e()  =  6.105  exp  (  25.22  — - —  -  5.31  In  ^ 


273.2 


4.  For  thermally  neutral  and  stable  conditions  (0  <  Rib  <  1),  the  evaporation 
duct  height,  6,  is  calculated  by 

(a)  6  =  0  for  AN  >  0 


or 


(b)  6  = 


A  Nn 


b\  B  -  ANp  a/L’ 

or,  if  the  result  of  (b)  is  such  that  6  <  0  or  6/  L’  >  1 ,  then 

,  ,  s  _  AjVp(I  +  a)  -  b,  a  z, 
b i  In  (z,/z0) 

For  stable  conditions,  the  function  B  is 


B  =  In  I  —  )  + 


a  z  i 


L' 


5.  For  thermally  unstable  conditions  (Rih  <  0),  evaporation  duct  height,  6,  is 
calculated  by 


6  = 


b  | 

AAL 


4/3  /*,  2?\3 


L'  \  AN, 


and  the  function,  B,  is 


B  -  In  (  ~  \  -  th 

zo 
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The  universal  function,  t A,  is  determined  by 


-e- 

II 

1 

4^ 

L/i 

Z  1 

>  -0.01 

i: 

L' 

t/r  =  I0[  ' 

02  log  ( 

r,  /.  )  -  0.69  ] 

-0.01 

V 

1^ 

IV 

L' 

II 

o 

.776  log i 

(  /.')  +  0.t06  ] 

-0.026 

V 

l-i* 

IV 

L 

lb  =  io[0 

,630  log  | 

:  -'|  /.’)  *0  )6  J 

-0.1 

L' 

o 

"o 

II 

S3- 

.414  log  l 

:  n  *  o.i6  ] 

-1 

V 

I-4 

IV 

L 

=  2 

-i 

<  -2.2 

L' 

-0.026 


0.1 

-1 

-2.2 


6.  The  maximum  value  of  evaporation  duct  height  is  limited  to  40  meters. 
Remaining  variables  and  constants  are  defined  as  follows: 


ANp  =  Na  -  Ws  ,  potential  refractivity  difference  between  the  air  and  the  sea 
Zi  =  6  meters,  reference  altitude 

b  |  =  -0. 125  m  1  ,  critical  gradient  of  potential  refractivity 
z0  =  1.5  x  10  4  m  ,  surface  roughness  parameter 
fi  -  4.5  ,  coefficient  in  the  unstable  function 
a  -  5.2  ,  coefficient  in  the  stable  function 


L'  - 


IQ--,  r, 

Rib 


m  ,  Monin-Obukhov  length,  where 


Ve  =  0.05 

T,  =  0.065  +  0.004  Rib 
r,  =  0.109  +  0.367  Rib 
re  =  0.155  +  0.021  Rih 


Rih  <  -3.75 
3.75  <  Rib  <  -0.12 
-0.12  <  Rib  <0.14 
0.14  <  Rib 


The  dominant  factor  in  determining  evaporation  duct  height  is  the  difference  in 
potential  refractivity  between  the  air  at  the  reference  altitude  and  the  air  at  the  sea  surface. 
Errors  in  air  temperature  measurement  caused  by  conductive  and  radiative  heating  effects 
have  been  shown  to  strongly  affect  the  duct  height  calculation.12  To  minimize  the  sensitivity 
of  the  evaporation  duct  algorithm,  an  additional  test  is  applied  whenever  Ta  -  Ts  >  -  1. 
Evaporation  duct  height  is  calculated  for  Ta  -  Ts,  ( 50 ),  and  for  Ta  -  T,  -  1,  (5  (),  with  Ts, 
u,  and  RH  unchanged.  Then,  if  <50  >  8  , ,  the  value  of  8  t  is  the  evaporation  duct  height; 
otherwise,  the  evaporation  duct  height  is  calculated  as  defined  in  steps  1  through  6  above. 
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4.1.2  Measurements  and  Procedures 


The  bulk  meteorological  measurements  required  for  the  evaporation  duct  calcula¬ 
tion  are  wind  speed,  air  and  sea  temperature,  relative  humidity,  and  altitude  of  the  measure¬ 
ments.  Since  gradients  in  the  surface  layer  generally  become  small  at  altitudes  more  than  a 
few  meters,  the  influence  of  measurement  altitude  is  neglected  as  long  as  measurements  are 
made  at  least  6  meters  above  the  surface.  Additionally,  the  influence  of  sea  waves  on  the 
measurements  is  usually  restricted  to  altitudes  less  than  three  times  the  wave  height  (Ref.  7, 
p.  123,  and  Ref.  13),  and  the  altitude  of  6  meters  or  greater  minimizes  most  wave  effects. 

4. 1.2.1  Wind  Speed 

Wind  speed  can  be  measured  with  the  typical  cup  or  aerovane  anemometer.  The 
anemometer  should  be  mounted  so  as  to  minimize  the  disturbance  to  the  air  flow  caused  by 
the  observing  platform.  On  a  ship  or  a  boat,  the  anemometer  should  be  mounted  as  far 
forward  and  as  high  as  possible;  the  measured  relative  wind  must  be  corrected  for  ship’s 
course  and  speed  to  determine  true  wind.  On  piers  or  other  coastal  structures,  the  anemo¬ 
meter  should  be  mounted  so  it  faces  into  the  prevailing  wind  and  is  at  least  2  meters  above 
the  top  of  the  structure.  Hand-held  anemometers  can  also  be  used  successfully  if  carried  to 
the  windward  side  of  the  platform. 

The  wind  speed  should  be  averaged  over  a  period  of  5  to  10  minutes.  This  will  aver¬ 
age  out  the  turbulent  fluctuations  in  wind  speed.  Figure  4  is  an  overwater  power  spectrum 
of  wind;  a  frequency  of  2  *  10  3  to  3  x  10  3  Hz  represents  an  averaging  period  over  approx¬ 
imately  5  to  10  minutes  and  avoids  the  turbulent  peak  at  0.01  Hz. 14  15 


Figure  4.  Byshev  and  Ivanov14  composite  average  wind  velocity  power 
spectrum  for  the  marine  atmospheric  surface  layer  from  data  taken  over 
the  South  Atlantic.  (From  Ref.  15.) 
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4. 1.2. 2  Air  Temperature  and  Relative  Humidity 

Measurement  of  air  temperature  and  relative  humidity  can  be  made  at  any  well- 
exposed  site  above  an  altitude  of  6  meters  with  reasonable  precautions  taken  to  minimize 
platform-induced  conductive  and  radiative  heating.  An  aspirated  psychrometer  with  wet 
and  dry  bulb  thermometers  is  a  convenient  portable  instrument.  Thermometers  or  other 
types  of  temperature  and  humidity  probes  can  be  installed  in  a  thermal  screen  mounted  in  a 
location  exposed  to  the  ambient  air.  Temperature  and  moisture  have  spectral  properties 
similar  to  those  shown  in  Fig.  4. 16  However,  since  temperature  and  humidity  gradients 
above  6  meters  in  the  typical  overwater  surface  layer  are  quite  small,  averaging  times  for 
these  parameters  can  be  3  to  5  minutes. 

4.1.2.3  Sea  Surface  Temperature 

The  simplest  measurement  of  sea  surface  temperature  is  made  with  a  small  rubber 
bucket  (or  any  container  with  a  small  thermal  mass)  and  an  accurate  thermometer.  More 
durable,  specifically  designed  water  thermometer  buckets  are  available  which  provide  a 
frame  for  a  thermometer  and  a  small  reservoir  (bucket)  to  retain  water  around  the  ther¬ 
mometer  bulb.  These  devices  are  used  to  sample  the  water  over  a  finite  depth  (~0.3  meter). 
Under  most  conditions,  the  wind  is  sufficient  to  mix  the  upper  levels  of  the  water  column 
(>2  meters).  However,  on  cloudless  days  with  strong  solar  heating  and  light  winds,  the 
surface  temperature  can  vary  considerably  from  the  underlying  water.  In  these  instances, 
care  must  be  taken  to  get  a  representative  sample  in  the  upper  0.3  meter  of  the  sea. 

Infrared  radiometer  sensors  provide  a  remote  means  to  sense  the  sea  surface 
temperature.  Although  considerably  more  expensive  than  the  buckets,  they  do  have  the 
advantage  of  sensing  only  the  sea  surface.  If  this  sort  of  instrument  is  used,  it  should  have  a 
means  to  compensate  for  reflected  sky  radiation. 

Further  information  on  instruments  and  procedures  can  be  found  in  Ref.  7,  15,  and 
17  19.  Table  3  shows  the  minimum  specifications  for  surface  measurements. 

Table  3.  Minimum  specifications  for  surface  measurements. 


Parameter 

Precision 

Range 

Wind  direction 

±5  deg 

5  to  360  deg 

Wind  speed 

±2  kt 

0  to  30  kt 

Air  temperature 

±0.5°  C 

-10°  to  40° C 

Relative  humidity 

±59 

0  to  100% 

Sea  temperature 

±0.5°  C 

-1°  to  40° C 

Wave  height 

±0.5  m 

0  to  6  m 

4.2  OCEAN  WAVES 

The  state  of  the  sea  affects  the  sea  clutter  return  for  radar  systems.  The  disturbance 
of  the  sea  surface  varies  in  scale  from  ripples  through  short-crested  (wind)  waves  to  long- 
wavelength  swells.  Ripples  (capillary  waves)  are  the  immediate  response  of  the  sea  surface 
to  the  wind  stress.  Short-crested  waves  develop  as  the  average  wind  continues  to  exert  a 
stress  on  the  sea  surface.  If  the  wind  acts  over  a  large  enough  area  (fetch)  for  a  sufficient 
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period  of  time  (duration),  a  steady-state  condition  is  reached  (fully-arisen  sea).  Swells  are 
long-wavelength  periodic  waves  formed  by  the  dispersion  of  wind  waves  as  they  propagate 
away  from  the  area  of  generation.  The  local  sea  state  is  determined  by  the  superposition  of 
the  locally  generated  ripples  and  wind  waves  and  the  swell  propagating  through  from 
distant  (possibly  multiple)  sources.  Wind-wave  height  is  well-correlated  to  the  average  wind 
speed  in  the  open  ocean.  Nearby  land  masses  or  major  ocean  currents  (e.  g..  Gulf  Stream) 
may  affect  this  correlation  by  creating  lower  or  higher  wave  heights. 

To  quantify  the  sea  conditions  for  sea-clutter  modeling,  it  is  necessary  to  measure 
wind  speed,  as  discussed  above,  wind  direction,  and  significant  wave  height.  Significant 
wave  height  is  defined  as  the  average  of  the  highest  one-third  of  the  waves  and  is  correlated 
to  what  a  careful  observer  would  visually  judge  as  the  predominant  wave  height.20-21  The 
best  way  to  estimate  the  significant  wave  height  is  by  observing  the  variation  of  the  height 
of  the  sea  surface  against  a  ship's  freeboard  or  a  coastal  structure’s  piling.  More  precise 
estimates  can  be  made  by  analyzing  wave  recorder  records  or  buoy  data  if  available. 

Wind  direction  can  be  measured  by  a  wind  vane  as  part  of  a  wind  measuring  system 
or  estimated  by  an  observer.  Wind  direction,  like  wind  speed,  should  be  averaged  over  5  to 
10  minutes.  Relative  wind  direction  determined  aboard  a  moving  ship  must  be  corrected 
for  ship’s  course  and  speed. 

Table  3  includes  minimum  specifications  for  ocean  wave  measurements. 
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5.0  EQUIPMENT  REFERENCES 


The  following  sample  listing  of  some  manufacturers  and  retailers  is  provided  as 
possible  sources  for  various  types  of  instruments  and  equipment.  This  listing  is  not  all- 
inclusive  nor  is  it  intended  to  be  an  endorsement  of  those  companies  or  products. 

•  Atmospheric  Instrumentation  Research  Inc.,  8401  Baseline  Road  A,  Boulder, 
CO  80303:  USA:  radiosonde  systems,  aerostats. 

•  Belfort  Instrument  Company,  727  South  Wolfe  Street,  Baltimore,  MD  2 1 23 1 : 
USA:  meteorological  instruments. 

•  Everest  Interscience  Inc.,  1 120  South  Raymond  Avenue,  Fullerton,  CA  92631: 
USA:  infrared  thermometers. 

•  Kahlsico  International  Corporation,  P.O.  Box  947,  El  Cajon,  CA  92022:  USA: 
meteorological  and  oceanographic  equipment;  water  bucket  thermometer. 

•  Qualimetrics  Inc.,  P.O.  Box  41039,  Sacramento,  CA  95841,  Representatives 
Science  Associates,  P.O.  Box  230,  Princeton,  NJ  08542,  and  Qualimetrics  LTD, 
%  Gold  Crest,  Prospect  Business  Center,  Dundee  Technology  Park,  Unit  No. 

24,  Dundee,  Scotland  DD21TY,  U.K.:  radiosonde  systems  and  meteorological 
instruments. 

•  Rotronic  Instrument  Corporation,  7  High  Street,  Huntington,  NY  1 1743:  USA: 
meteorological  probes. 

•  Vaisala  Oy,  PL  26,  SF-00421,  Helsinki  42,  Finland,  with  subsidiaries  Vaisala 
Inc.,  2  Tower  Office  Park,  Woburn,  MA  01801  and  Vaisala  (UK)  LTD,  1 1 
Billing  Road,  Northampton  NN1  5AW,  U.K.:  radiosonde  systems  and  surface 
meteorological  instruments. 

•  VIZ  Manufacturing  Company,  335  East  Price  Street,  Philadelphia,  PA  19144- 
5782  USA:  radiosonde  systems. 
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6.0  GLOSSARY 


The  semantics  of  the  various  terms  that  are  used  in  describing  the  performance  of 
instruments  necessitates  defining  the  terms  as  they  are  used  in  this  document.  These  defini¬ 
tions  follow  Brock.18 

accuracy  —  used  in  this  document  only  as  a  general,  descriptive  term  of  an  instru¬ 
ment’s  performance. 

lag  —  the  time  difference  between  when  a  specific  value  of  output  is  obtained  and 
when  that  same  value  was  previously  obtained  by  the  input. 

precision  —  a  statistical  description  of  the  amount  of  deviation  from  the  most  likely 
value  for  the  measured  quantity. 

range  —  the  interval  of  the  measured  quantity  over  which  the  instrument  is 
designed  to  be  sensitive. 

resolution  —  the  smallest  variation  in  the  input  that  causes  a  detectable  change  in 
the  output. 
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Appendix  A 


DEVELOPMENT  OF  THE  EVAPORATION  DUCT  HEIGHT  FORMULATION 


References  A-l  to  A-3  provide  a  detailed  development  of  the  evaporation  duct 
height  formulation  by  Jeske.  A  general  development  of  Jeske’s  formulation  is  included  here 
along  with  assumptions  and  rationale  for  the  implementation  in  the  computer  program  list¬ 
ing  in  Appendix  D. 

Jeske  defined  the  gradient 


di  SL 

—  =  - ‘ -  d> 

dz  p  ku,  (z  +  r0) 


(A-l) 


where  i  is  a  scalar  property  of  the  atmosphere,  St,  is  the  vertical  flux  of  i,  z  is  height  in 
meters,  k  is  von  Karmen’s  constant  (0.4),  u.  is  the  friction  velocity  in  m/s,  z0  is  the  rough¬ 
ness  parameter  in  meters,  and  d>  is  a  stability-dependent  function  (equal  to  1  at  thermally 
neutral  conditions,  i.  e.,  air  temperature  equal  to  sea  surface  temperature).  For  stable 
conditions  (air  warmer  than  water),  the  stability  function  was  taken  to  be  the  logarithmic- 
linear  model  proposed  by  Monin  and  Obukhov: 


<t> 


-  I  +  a 


L' 


( A-2) 


where  a  is  taken  to  be  5.2  and  L’  is  the  gradient  form  of  the  Monin-Obukhov  scaling 
length  corrected  for  stability.  For  unstable  conditions  (water  warmer  than  air),  the  KEYPS 
relationship 


d>4 


(A-3) 


with  a  -  4.5  is  used.  To  evaluate  z/ L\  the  bulk  Richardson  number  is  used  as 


g±Tz 
Tu2  r 


(A-4) 


where  g  -  9.8  m/s2,  c  is  height  in  meters.  7"  is  air  temperature  in  Kelvin,  AT  is  air  tempera¬ 
ture  minus  sea  temperature,  u  is  wind  speed  in  m/s,  and  T  is  the  profile  coefficient.  Rih 
involves  approximating  gradients  with  differences  and  assumes  7"*  8  in  the  surface  layer. 
The  profile  coefficient  for  neutral  conditions  is 


I 


and  is  approximately  0.1  for  typical  observation  heights,  z{,  and  surface  roughness  parame¬ 
ter  values,  z0,  of  0.00015  meter.  Similarly,  L'  is 


Tu2  re 

L  -  - - 

gM 


(A -6) 


A-l 


or.  in  terms  of  Rih 


where  T,  is  an  empirical  profile  coefficient  for  nonneutral  conditions  (Fig.  A-l).  r?  is  eval¬ 
uated  by  a  straight-line-segment  fit  to  the  data  in  Fig.  A-l. 


Figure  A-1.  The  profile  coefficient,  as  a  function  of  bulk  Richardson 
number,  Rib.  (From  Jeske.A  l  10  A'3) 


Returning  to  Eq.A-1,  assume  the  flux,  Sr  is  constant  with  height  and  integrate 
from  the  sea  surface  to  z(: 


1 1  -  t0  -  Ai  |  -  j  -  dz 


P  K  U,  0  Z  +  Zf 


Thus,  the  flux  term  is 


P  x  ut  z,  0 

J  -  dz 

0  2  + 

and  substituting  in  Eq.  A-l  yields 


(A-8) 


(A-9) 


di  Ai  (p 
dz  Biz  +  z0) 

where 


(A- 10) 


A-2 


C  Z 1  <t> 

5=/  -  dz  ( A-l  1) 

0  2  +  ^0 

Taking  the  scalar  i  as  potential  refractivity,  )V  the  height,  6,  at  which  the  critical  gradient 
of  potential  refractivity,  6,  =  -0.125  m  ',  occurs  is  found  from  Eq.  A-10  as 


8  = 


btB-  AN'p  a/L' 


where  d>  for  stable  conditions  is  determined  from  Eq.  A-2  and  B  is  evaluated  as 


( A- 1 2) 


In  —  + 
:o 


( A- 1 3) 


Jeske  recommends  that  6/ L'  should  not  exceed  I.  If  this  occurs,  then  duct  height  should  be 
recalculated  with  8/ L’  -  I . 

For  unstable  conditions,  solving  Eq.  A-10  for  <t>  in  terms  of  Np,  substituting  in  Eq. 
A-3  with  z  =  6,  and  solving  yields 


5  = 


/  b,  B 

\  AN p 


4  a  /  b  |  B  \  ?  1/4 

\  AN p  I 


where  B  is  evaluated  as 


B=  In  ~  ) 

‘0  L 

4i(z\  I  L')  has  the  analytic  formA'4 


i/r  =  1  -  <f>  -  3  lnd>  +  2  In 


7T 


+  2  tan  1  4> - +  In 

2 


(A-14) 


( A- 1 5) 


(A- 16) 


Use  of  Eq.  A-16  requires  the  iterative  solution  of  Eq.  A-3  for  </>  with  a  given  z/Z.'.  t jj  can  be 
more  quickly  determined  by  a  straight-line-segment  fit  to  the  K.EYPS  profile  of  Fig.  A-2 
adapted  from  Lumley  and  Panofsky.A_5 
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CRITICAL  GRADIENT  OF  POTENTIAL  REFR ACTIVITY 

The  equation  for  refractivity  in  terms  of  pressure,  P.  in  millibars,  temperature,  T \  in 
Kelvin,  and  water  vapor  pressure,  e,  in  millibars  is 

77.6  P  <  e 

N  =  -  +  3.73  x  105  —  (B-l) 

T  T2 

For  application  in  the  atmospheric  surface  layer,  potential  refractivity 
77.6  P„  .  en 

N„  =  - -  +  3.73  x  10s  -4  (B-2) 

p  8  82 


is  the  conservative  property  where  PQ  is  a  reference  pressure  level  (taken  to  be  1000  mb),  8 
is  potential  temperature,  and  ep  is  potential  water  vapor  p-sssurc. 

From  geometric  optics,  the  critical  gradient  required  for  trapping  is  that  which 
yields  a  ray  curvature  equal  to  the  earth’s  curvature 

dN  106 

-  =  _  -  =  -0.157  N/m  (B-3) 

d:  a 


where  a  is  the  earth’s  radius,  6371  *  103  m.  To  relate  the  gradient  of  potential  refractivity  to 
the  gradient  of  refractivity,  take  the  total  derivatives  of  Np  (8,  ep )  and  N(P,  T,  e)  with 
respect  to  z.  Assuming  8  *=  T  and  ep  e  in  the  surface  layer  yields 


dNp  _  dN  dN  dP 
dz  dz  dP  dz 


(B-4) 


The  partial  derivative  dN/dP  =  77.6/  T  varies  from  0 .'*8  to  0.26  mb  1  over  the  temperature 
range  of  0°  to  30° C.  The  derivative  dPjdz  can  be  evaluated  by  making  the  hydrostatic 
approximation 


dP 

dz 


=  -P  g~  - 


Pg 

R  T 


(B-5) 


where  p  is  density,  g  is  acceleration  of  gravity  (9.8  m/s2),  P  is  pressure  (1000  mb),  R  is  the 
individual  gas  constant  for  dry  air  (2.87  *  106  erg/g/°K),  and  7" is  temperature  (°K).  This 
yields  a  variation  of  dPjdz  from  -0.12  to  -O.f  1  mb/m  over  the  temperature  range  of  0°  to 
30°C.  Taking  standard  temperature  (15°C),  the  critical  gradient  for  trapping,  in  terms  of 
potential  refractivity,  is 

dN 

—P--  -  0.157  -  (0.27)(-0. 1 2)  =  -  0.125  m  1  (B-6) 

dz 
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IMPLEMENTING  A  MODIFIED  EVAPORATION 
DUCT  HEIGHT  CALCULATION 

Reference  C-l  showed  that  a  bias  in  air-sea  temperature  difference  toward  thermal 
stability  resuited  in  spuriously  high  calculated  evaporation  duct  heights.  A  modification  to 
the  duct  height  calculation  that  detected  and  compensated  for  the  air-sea  temperature 
difference  bias  was  proposed  based  on  the  variation  of  duct  height  as  a  function  of  air-sea 
temperature  difference  at  a  constant  relative  humidity.  To  implement  the  algorithm,  the 
value  of  either  the  air  temperature  or  sea  temperature  may  have  to  be  changed  to  yield  a 
modified  air-sea  temperature  difference.  A  comparison  of  climatological  air  and  sea 
temperatures  from  Marsden  square  data0'2  and  ocean  data  buoys0'3  showed  better  agree¬ 
ment  between  sea  temperature  means  than  between  air  temperature  means.  Air  temperature 
means  from  the  Marsden  square  data  also  indicate  a  diurnal  variation  (~  1°C)  that  the 
buoy  data  do  not  show.  Mean  air  temperatures  for  3-hour  intervals  for  data  buoys  show 
diurnal  variations  of  only  a  few  tenths  of  a  degree.*  This  gives  more  weight  to  holding  sea 
temperature  constant  and  varying  air  temperature  to  obtain  the  appropriate  air-sea 
temperature  difference. 

The  next  question  that  arises,  then,  is  whether  or  not  it  is  reasonable  to  hold  rela¬ 
tive  humidity  constant  while  changing  air  temperature.  Reed  reported  finding  systematic 
errors  in  air  temperature  measurements  aboard  a  ship  but  did  not  find  systematic  errors  in 
relative  humidity  measurements.04  Experience  with  psychrometric  measurements  in  a 
thermal  screen  during  the  Ku-Band  Surface  Surveillance  Project0'5  sited  on  Point  Loma  in 
San  Diego  tends  to  support  Reed’s  findings.  Recommendations  of  a  working  group  on  the 
computation  of  global  air-sea  flux  climatology  to  use  both  daytime  and  nightrme  humidi¬ 
ties  but  only  nighttime  air  and  sea  temperature  lend  support  to  the  assumption  of  negligible 
bias  in  the  relative  humidity  measurements.0"6 

The  implementation  of  the  modified  evaporation  duct  height  calculation  assumes 
that  air  temperature  can  be  varied  to  obtain  the  appropriate  air-sea  temperature  differences 
while  sea  temperature,  wind  speed,  and  relative  humidity  are  held  constant.  Figure  C-l  is  a 
flow  chart  of  the  process. 


*D.B.  Gilhousen,  National  Oceanic  and  Atmospheric  Administration  Data  Buoy  Center,  personal 
communication,  January  1984. 
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Figure  C-1 .  Flow  chart  for  the  modified  evaporation  duct 
calculation  based  on  air- sea  temperature  difference  (ASTD). 
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FORTRAN  PROGRAM  TO  CALCULATE  EVAPORATION  DUCT  HEIGHT 


c 

c  name:  delta. f 

c  routine  controls  the  execution  of  the  'getdelta' 

c  function  for  testing  the  calculations  of  the  evaporation 

c  duct  height  (delta) . 

c  language:  fortran  iv  ('77) 

c  machine:  hp  520  (unix) 

c  call  sequence: 

c  MAIN  routine 


c 

glossary : 

c 

airt : 

air  temperature,  deg  C 

c 

seat : 

sea  temperature,  deg  C 

c 

rh: 

relative  humidity,  % 

c 

u : 

wind  speed,  knots 

c 

subroutines : 

c 

getdelta 

( function) 

c 

evap 

c  rev  date  description 

c  0.00  101987 

c  ****************************************************** 

c  prompt  the  user  for  entries,  call  'getdelta'  to 

c  compute  the  evaporation  duct  height 

c  ****************************************************** 

write(*,9000) 

9000  format(  "Compute  the  evaporation  duct  height  from"/ 


* 

"entries  of:"/ 

* 

"  air  temperature,  deg  C"/ 

* 

"  sea  temperature,  deg 

C"/ 

* 

"  relative  humidity,  % 

and"/ 

★ 

"  wind  speed,  knots"/) 

write(7 , 9010) 

9010  format(" . . air ... | . sea ... | . . rh ... | ... u. . . | . delta . | "/ 
*  "  temp  j  temp  j  %  j  knots  j  meters") 

do  while ( . true . ) 

write (*,'( "enter  air  temp  (C) :  ",$)') 
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read(*,*)  airt 


write(*, ' ("enter  sea  temp  (C) :  ",$)') 
read(*,*)  seat 

write(*. ' ("enter  relative  humidity  (%):  ",$)') 
read(*,*)  rh 

write(*, ' ("enter  wind  speed  (knots):  ",$)') 
read(*,*)  u 

delta-get_delta(airt ,  seat,  rh,  u) 
write(7 , 9010) 

write (7 , 9020)  airt, seat, rh.u, delta 
9020  format(5f8 . 1) 

enddo 

end 
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c 

c  ****************************************************** 

c  name:  getdelta.f 

c  routine  calculates  the  evaporation  duct  height  using 

c  the  paulus  formulation  "Practical  Application  of  the 

c  IREPS  Evaporation  Duct  Model",  NOSC  TR  966,  Jun  84 

c  language:  fortran  iv  ('77) 

c  machine:  hp  520  (unix) 

c  call  sequence: 

c  delta-get_delta(airt ,  seat,  rh,  u) 

c  glossary: 

c  airt 

c  seat 

c  rh : 

c  u : 

c  subroutines: 

c  evap 

c  rev  date  description 

c  0.00  101987  modify  'duct63'  program  code  to 

c  make  a  stand-alone  utility. 

c  ***********  ****************************************** 

function  get_delta(airt , seat , rh,u) 

c  ****************************************************** 


c  check  wind  speed,  convert  to  deg  K,  and  compute 

c  vapor  pressures 

if  (u  .gt.  .01)  then 

wvel=»u 

dtok-273.2 

tak-airt+dtok 

tsk-seat+dtok 

esw-6 . 105*exp(25 . 22*seat/tsk- 5 . 31*alog( tsk/dtok) ) 
es  -6 . 105*exp(25 . 22*airt/tak- 5 . 31*alog( tak/dtok) ) 
e  -rh*es/100. 

c  calculate  delta 

tl-tsk-1.0  [modified  air  temp 


air  temperature,  deg  C 
sea  temperature,  deg  C 
relative  humidity,  % 
wind  speed,  knots 
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if(tak  .gt.  tl)  then 
el-rh*esw/100 . 

call  evap(wvel , tsk , tsk , el , esw, delO) 

esl-6 . 105*exp(25 .22*(seat-1.0)/tl-5. 31*alog( tl/dtok) ) 
el-rh*esl/100 . 

call  evap(wvel , tl , tsk , el , esw, delta) 
if (delO-delta  .It.  0.)  then 

call  evap(wvel , tak, tsk, e , esw, delta) 
endif 
else 

call  evap(wvel , tak, tsk , e , esw, delta) 
endif 


else 


delta-0 . 
end  if 

get_delta-aminl (delta , 40 . ) 

return 

end 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


***************************************************** 


Name :  EVAP 


Routine  calculates  the  evaporation  duct  height 
(in  metres)  following  the  method  of  Jeske  "The  State 
of  Radar  Range  Prediction  over  Sea",  Tropospheric 
Radio  Wave  Propagation  -  Part  II,  NATO-AGARD,  Feb  71 

Language:  FORTRAN  IV  ('77) 

Machine  :  Univac  1100/82 


CALL  Evap(wvel , tak, tsk, e , esw, delta) 


wvel : 
tak: 
tsk  : 
e  : 
esw  : 
delta 


Wind  speed  knots 
Air  temperature  in  Kelvin 
Sea  temperature  in  Kelvin 
Vapor  pressure  mbs 
Surface  vapor  pressure  mbs 
Duct  height  m  RETURNED 


Rev  Date  Description 

0.0  052984 


***************************************************** 


subroutine  evap(wvel , tak , tsk, e , esw , del ta) 
real  lnzlzO 


delta-0 . 

z0-l . 5e-4 

zl-6 . 

beta-4 . 5 

alpha-5 . 2 

Inzlz0-alog(zl/z0) 

bl--0. 125 


! default  duct  height 
!  roughness  height 
!  reference  height  (m) 

!  constant 
!  constant 
!  useful 

(potential  refractivity  trapping  grad. 


rib-369 . *zl*(tak- tsk)/( tak*wvel*wvel) 
if(rib  .gt.  1.)  rib-1. 

compute  potential  refractivity  difference 

delna-77 . 6/ tak* (1000 . +4810 . *e/tak) 
delnO— 77 . 6/tsk*( 1000 . +4810 . *esw/tsk) 
de ln-de lna -delnO 
if (dein  . ge .  0.)  return 


compute  the  gamma  function 
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c 


if (rib  .le.  -3.75)  then 
gamma-0 . 05 

else  if(rib  .le.  -0.12)  then 
gamma-0 . 065+rib*0 . 004 
else  if(rib  . le .  0.14)  then 

gamma-0 . 109+rib*0 . 367 

else 

gamma-0 . 155+rib*0 . 021 
end  if 
c 

olp-rib/(10.*zl*gamma)  ll./L' 

zlolp-zl*olp 
c 

c  . 

c  compute  Psi  function  for  unstable  conditions 

c 

if (rib  .It.  0.)  then 

if(zlolp  .ge.  -0.010)  then 
psi--4 . 5*zlolp 

else  if(zlolp  . ge .  -0.026)  then 

psi-10 . **(1 . 020*alogl0( -zlolp) +0.690) 
else  if'zlolp  .ge.  -0.100)  then 

psi-10 . **(0 . 776*alogl0( - zlolp) +0 . 306) 
else  if(zlolp  .ge.  -1.000)  then 

psi-10  **(0 . 630*alogl0( -zlolp )+0. 160) 
else  if (zlolp  .ge.  -2.200)  then 

psi-10.**(0.414*alogl0(-zlolp)+0. 160) 

else 

psi-2 . 000 
end  if 
c 

c  Compute  unstable  delta 

c 

b-lnzlzO-psi 

btemp-b*bl/deln 

d-btemp**4  -  4 . *beta*olp*btemp**3 
if(d  .gt.  0.)  delta-d**( -0 . 25) 
c 

c  Stable  and  neutral  conditions 

c 

else 

b-lnzlz0+zlolp*alpha 

delta-deln/(b*bl  -  alpha*deln*olp) 

if(delta  .It.  0.  .or.  delta*olp  .gt.  1.)  then 

delta-(deln*( 1 . +alpha) -bl*alpha*zl)/(bl*lnzlzO) 
end  if 
end  if 
c 

c  Return  to  caller 

c 

return 

c 

end 
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TEST  CASES  FOR  EVAPORATION  DUCT  HEIGHT  CALCULATION 

For  the  specified  inputs  of  air  temperature  (°C),  sea  temperature  (°C),  relative 
humidity  (%),  and  wind  speed  (knots),  the  evaporation  duct  height  (6)  should  be  calculated 
correctly  to  within  ±  0.1  meter. 


Air  Sea  RH  u  8 

Temp(°C)  Temp(°C)  (%)  (knots)  (m) 


0.0 

15.0 

10. 

1. 

3.9 

0.0 

15.0 

10. 

20. 

16.8 

0.0 

15.0 

75. 

1. 

2.8 

0.0 

15.0 

75. 

20. 

12.3 

0.0 

15.0 

100. 

1. 

2.4 

0.0 

15.0 

100. 

20. 

10.4 

0.0 

16.0 

10. 

1. 

4.0 

0.0 

16.0 

10. 

20. 

17.3 

0.0 

16.0 

75. 

I. 

3.0 

0.0 

16.0 

75. 

20. 

13.0 

0.0 

16.0 

100. 

1. 

2.6 

0.0 

16.0 

100. 

20. 

11.2 

15.0 

15.0 

10. 

1. 

8.7 

15.0 

15.0 

10. 

20. 

37.0 

15.0 

15.0 

75. 

1. 

3.6 

15.0 

15.0 

75. 

20. 

13.9 

15.0 

15.0 

100. 

1. 

0.0 

15.0 

15.0 

100. 

20. 

0.0 

15.0 

16.0 

10. 

1. 

9.1 

15.0 

16.0 

10. 

20. 

38.8 

15.0 

16.0 

75. 

0 

0.0 

15.0 

16.0 

75. 

1. 

3.8 

15.0 

16.0 

75. 

20. 

14.6 

15.0 

16.0 

100. 

1. 

0.9 

15.0 

16.0 

100. 

20. 

2.6 

30.0 

15.0 

10. 

1. 

8.7 

30.0 

15.0 

10. 

20. 

37.0 

30.0 

15.0 

75. 

1. 

3.6 

30.0 

15.0 

75. 

20. 

13.9 

30.0 

15.0 

100. 

1. 

0.0 

(Cont’d) 
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Air 

Temp  (°C) 

Sea 

Temp  (°C) 

RH 

(%) 

u 

(knots) 

6 

(m) 

30.0 

15.0 

100. 

20. 

0.0 

30.0 

16.0 

10. 

1. 

9.1 

30.0 

16.0 

10. 

20. 

38.8 

30.0 

16.0 

75. 

1. 

3.8 

30.0 

16.0 

75. 

20. 

14.6 

30.0 

16.0 

100. 

1. 

0.0 

30.0 

16.0 

100. 

20. 

0.0 
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M-PROFILES  IN  THE  SURFACE  LAYER 

Propagation  models  require  a  refractivity  profile  to  make  their  calculations.  The 
procedures  specified  in  the  main  text  yield  only  a  scalar  duct  height.  In  a  manner  similar  to 
Jeske,*  an  equation  for  M  versus  height  can  be  developed.  The  equation  for  potential 
refractivity,  Np ,  is 

77.6  P„  .  en 

S  =  - "  +  3.73  «  105  -4  (F-l) 

p  d  0 2 


where  Pa  is  a  reference  pressure  (1000  mb),  6  is  potential  temperature  (°K),  and  ep  is 
potential  water  vapor  pressure  (mb).  The  total  derivative  of  N  with  respect  to  altitude. 


is 


dNp  _  dNp  d6  dNp  dep 
dz  dd  dz  de„  dz 


(F-2) 


The  partial  derivatives  in  Eq.  F-2  are  evaluated  as 


dNp 

dd 


77-6  P0 
02 


-  7.46  *  105 


dNp  _  3.73  *  105 
dep  e2 


(F-3) 


Over  the  range  of  0°  to  30°C,  dNp/dO  varies  from  -1.2  to  -1.6/°K;  dNp/den  varies  f,om 
5.0  to  4.1  mb  (for  potential  vapor  pressure  equal  to  80%  of  saturation).  Then,  for  15°C 
and  80%  relative  humidity.  Eq.  F-2  becomes 


dNp 

dz 


1.4 


dd  de 

—  +  4.5  _27 


dz  dz 


(F-4) 


The  potential  temperature  and  potential  water  vapor  profiles  in  the  surface  layer  are  found 
from 


dl~.  »(.!) 

dz  k(z  +  zfl)  \  I 7  / 


%« 


K(z  +  Zq) 


(F-5) 


♦Jeske,  H.,  “State  and  Limits  of  Prediction  Methods  of  Radar  Wave  Propagation  Conditions  over 
Sea,”  in  Modern  Topics  in  Microwave  Propagation  and  Air- Sea  Interaction,  A.  Zancla,  ed.,  Reidel 
Pub.,  1973. 
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where  0,  and  ep,  are  scaling  quantities,  k  is  von  Karmen’s  constant,  z  is  altitude,  z0  is  the 
aerodynamic  roughness  parameter,  L'  is  the  Monin-Obukhov  scaling  length  corrected  for 
stability,  and  4>{z/  L')  is  a  stability-dependent  function.  Equation  F-4  then  becomes 


dNp 

dz 


1.4 


6. 

k(z  +  z0) 


+  4.5 


k(z  +  z0) 


(F-6) 


The  scaling  quantities  for  a  given  set  of  conditions  are  determined  from  bulk  meteorologi¬ 
cal  measurements 


6,  -  K 


0(Zj)  -  0(zo) 


(F-7) 


V  =  * 


?p(z\)  -  ep(z o) 


where  z,  is  the  reference  altitude  and  ^(z,/L')  is  the  integrated  stability-dependent  func¬ 
tion.  Substitution  of  Eq.  F-7  into  Eq.  F-6  yields 


dN 


p 


dz 


(F-8) 


where  AiVp  is  the  bulk  air-sea  potential  refractivity  difference 

SNp  =  -  1 .4  [0(z, )  0(zo)]  +  4.5  \ep(zx)  -ep(z0)]  (F-9) 

The  evaporation  duct  height,  6,  is  the  height  at  which  the  gradient  of  potential  refractivity 
is  -0.125  m  ’.  Thus  Eq.  F-8  becomes 


0.125  = 


(5  +  z0) 


ANp 


(F-10) 


Solving  Eq.  F-10  for  ANp  and  substituting  into  Eq.  F-8  gives 


dNp 

dz 


-  -  0.125 


(6  +  ^0) 


*{t)  ,z*"> 


(F-ll) 


F-2 


But,  dNp/ dz  -  dM/d z  -  0.125,  and  the  A/-gradient  is 


dM  (8  +  z0)  *  \L’) 

-  =  0.125  -  0.125  - - —  - 

dz  (  8  \  (z  +  z0) 

*\T') 


(F-12) 


Integrating 


8  +  z0  fz  *\L') 

M(z)  =  M  +  0.125  z  -  0.125  - —  /  -  dz 

*(t)  ° 


(F-13) 


For  stable  conditions,  the  stability-dependent  function  has  the  form,  <t>  =  1  +  a  z/L’,  and 
Eq.  F-13  becomes 


o  /  z  +  zn  az 

M(z)  -  M  +  0.125  z  -  0.125  -  (In - -  +  — 

a8  \  z0  L' 


(F-14) 


where  z0  has  been  neglected  relative  to  8.  For  unstable  conditions,  the  stability-dependent 
function  has  the  form, 


d>4  -  4a  —  =  1 

U 


and  Eq.  F-13  becomes 


M(z)  -  Mp  +  0.125  z  -  0.125 


/  5 


z  +  z0 


(F-15) 


where  z0  has  again  been  neglected  relative  to  6.  d>  and  ^  can  be  determined  as  in  Appendix 
A.  It  has  been  found  through  numerous  propagation  experiments  and  model  simulations 
that  the  neutral  A/-profile  is  a  good  representation  of  practical  evaporation  duct  condi¬ 
tions.  For  neutral  conditions,  d>  =  1  and  ^  =  0,  and  Eq.  F-15  is  simply 


M(z)  -  M0  +  0.125  z  -  0.125  8  In 


z  +  z0 


(F-16) 


The  aerodynamic  roughness  parameter,  zn,  is  taken  to  be  constant  at  1.5  x  IO-4  meter. 
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